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ABSTRACT 

High-spatial resolution near-infrared (NIR) images of the central 24 x 24 arcsec 2 
(~2x 2kpc 2 ) of the elliptical galaxy NGC 1052 reveal a total of 25 compact sources ran- 
domly distributed in the region. Fifteen of them exhibit Ha luminosities an order of magnitude 
above the estimate for an evolved population of extreme horizontal branch stars. Their Ha 
equivalent widths and optical-to-NIR spectral energy distributions are consistent with them 
being young stellar clusters aged < 7 Myr. We consider this to be the first direct observation 
of spatially resolved star-forming regions in the central kiloparsecs of an elliptical galaxy. The 
sizes of these regions are < 11 pc and their median reddening is E(£> — V) ~ 1 mag. Accord- 
ing to previous works, NGC 1052 may have experienced a merger event about 1 Gyr ago. On 
the assumption that these clusters are spreaded with similar density over the whole galaxy, 
the fraction of galaxy mass (5 x 10~ 5 ) and rate of star formation (0.01 M Q yr _1 ) involved, 
suggest the merger event as the possible cause for the star formation we see today. 

Key words: galaxies: individual: NGC 1052 - galaxies: interactions - galaxies: nuclei - 
galaxies: star clusters - techniques: high angular resolution. 



1 INTRODUCTION 

About ~ 30% of the local early-type galaxies present signs of Re- 
cent Star Formation (RSF) in the last Gyr (Kaviraj et al. 2007). This 
result arises from the study of the integrated UV-optical colours of 
a large sample of early-type galaxies. A sizeable fraction of objects 
are revealed by their blue colours, NUV—r < 5.5 mag, despite of 
the fact that their optical colours place them in the optical "red se- 
quence" (u — r > 2.2 mag, Strateva et al. 2001). Given these blue 
NUV—r colours, a merger is likely the origin of the induced RSF 
(Kaviraj et al. 2009; Kaviraj 2010; Lopez-Sanjuan et al. 2010). 

In that scenario, one would expect the formation of new young 
cluster populations in an elliptical galaxy which presents blue UV- 
optical colours and a past merger event (Kaviraj 2010). However, 
most of the examples reported up to date exhibit intermediate- 
age clusters, between > 200 Myr and a few Gyr (e.g. NGC 1316, 
NGC 1380, NGC 1700, NGC 3610, NGC 5128 Carlson et al. 1999; 
Goudfrooij et al. 2001; Williams et al. 2010). In those cases, most 
of the age estimates are based on the comparison of cluster colours 
with stellar population models, and thus are strongly dependent 
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on the metallicity, extinction and IMF assumed in these models. 
Therefore, a scrutinized study of nearby ellipticals with blue UV- 
optical colours is particularly interesting to reveal the young cluster 
population at the earliest stage of their evolution, i.e. < 10 Myr. It 
is in this context that this letter is addressed. It presents the first 
detection of individually resolved young stellar clusters, < 7 Myr 
old, in a nearby elliptical galaxy, NGC 1052. To single out this pop- 
ulation, high spatial resolution imaging from UV to optical and in- 
frared (IR) are used. A complete spectral energy distribution (SED) 
is built for each of the clusters, their nature, ages and dust estimates 
have been addressed on the basis of their Ha emission. 

NGC 1052 is one of the nearest elliptical galaxies (E4, 
18 Mpc, 1" « 87 pc, lensen et al. 2003) with a stellar mass of 
Mi, « 10 11 M (derived from 2MASS 1 photometry, assuming 
M*/L K = 1.32M Q /L Q , Cole et al. 2001) and a low-level ac- 
tivity LINER nucleus (Heckman 1980). It presents two radio lobes 
(Cooper et al. 2007, Fig. 1) and a collimated ionized gas structure 
in Ha+[N II] (Pogge et al. 2000), both extended in the east-west 
direction, and a parsec-scale twin jet (PA-~60°) with an emis- 
sion gap between the brightest part located north-east (NE) and the 
counter jet in the south-west (SW) direction (Kadler et al. 2004). 



1 2 Micron All Sky Survey, see http : / /www . ipac . caltech . edu/ 
2mass/ 
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Figure 1. Left: SDSS g (blue), r (green) and i-band (red) colour composite image for NGC 1052. VLT/NaCo 24 X 24 arcaec 2 FOV is marked by a black 
square. Centre: Unsharp-masked VLT/NaCo Kv-band image, logarithmic contours for the natural Ks-band image (grey colour) and VLA contours at 1.4 GHz 
(black, beam ~ 2" X l'/5, Cooper et al. 2007). In blue, knots with EWHa ^ 50 A, in red those with < 50 A. The central white barred-dot marks the position 
of the nucleus and the direction of the parsec-scale twin jet. Right: SEDs for two "blue" (#1 and #19) and two "red" knots (#8 and #17). 



The radio lobes are also spatially coincident with the X-ray emis- 
sion region reported by Kadler et al. (2004). The distribution of the 
H I gas, with two tidal tails extended from NE to SW and also an 
infalling component, is interpreted as evidence of a merger event 
about 1 Gyr ago (van Gorkom et al. 1986, 1989). This scenario 
is further supported by the different rotation axes found for the gas 
and the stars (Plana & Boulesteix 1996) and by the presence of dust 
lanes along the N-S direction near the galaxy centre (Forbes et al. 
1990). Furthermore, NGC 1052 shows typical integrated colours of 
a "rejuvenated" early-type galaxy (following Kaviraj et al. 2007), 
having NUV—r — 4.9 mag, while the u — r — 2.8 mag colour 
places this galaxy in the optical red sequence (NUV from Ram- 
pazzo et al. 2007; u and r-band from SDSS"). 



2 OBSERVATIONS 

The brightness of the LINER nucleus permitted us to perform 
high-spatial resolution observations using the Very Large Tele- 
scope (VLT) and the Nasmyth Adaptive Optics System (NAOS) 
plus the Near-Infrared Imager and Spectrograph (CONICA). The 
VLT/NaCo dataset consist of near-infrared (NIR) J, Ks and L'~ 
band observations for the central 24 x 24arcsec 2 of NGC 1052 
(~2x2kpc 2 , Fig. 1), taken on 2005 December 3. The total 
exposure times are 900 s (/), 400 s (Ks) and 5.25 s (Z/), and 
the achieved FWHM resolution -measured in the most compact 
objects in these images- are <0'.'26, <0'.'12 and <0'.'12, re- 
spectively. Data reduction of the NaCo dataset was performed 
using the ECLIPSE package, provided by ESO, and includes 
sky subtraction, registration and combination of frames to ob- 
tain a final image per filter. The photometric calibration was 
based on standard stars taken along with the science frames. 

2 Sloan Digital Sky Survey, see http : / /www .sdss.org/ 



The relative photometric error is ~8%, ~ 11% and ~ 11% in 
the J, Ks and L' -bands, respectively. The VLT/NaCo dataset 
was complemented with the following Hubble Space Telescope 
(HST) imaging: WFPC/F555W (Wide Field Planetary Cam- 
era/V-band, A0.55p.rn AA0.05p.rn), WFPC2/F658N (including 
Hq+[Nii], A 0.659 pm AA 0.003 pjn), STIS/F28X50LP (Space 
Telescope Imaging Spectrograph, A 0.72 pm AA 0.11pm) and 
NICMOS/F160W (Near-Infrared Camera and Multi-Object Spec- 
trometer/rY-band, A 1.61pm AA0.12p.rn). The reduced images 
were retrieved directly from the scientific archive, and the photo- 
metric calibration was performed by using the PHOTFLAMB key- 
word value in the image header as a conversion factor. 

Image alignment of the whole set is based on the identifi- 
cation of two compact sources (knots) seen in all filters but in 
the L'-band. We avoided using the nucleus since dust obscura- 
tion may distort its position in the optical. For the L'-band, we as- 
sumed its position to match that in the Ks-band. We further used a 
1.4 GHz image with a beam resolution of ~ 2" x l'.'5 from VLA 
(Cooper et al. 2007). The signal-to-noise achieved in the NaCo J 
and fe-bands permitted us to distinguish a number of point-like 
sources in the FOV. In order to enhance the image contrast, dif- 
ferent techniques were applied, the best results being obtained by 
the unsharp-masking method (Sofue 1993). A median circular fil- 
ter with radius 0'.'12 was used to smooth the images, which were 
then subtracted from the original ones. The result for the Ks-band 
is shown in the background of Fig. 1 . The DAOPHOT finding algo- 
rithm (Stetson 1987) recovered, on the J and Ks unsharp-masked 
images, a total of 25 point-like sources over a 3<r threshold in the 
central 24 x 24arcsec 2 (circles in Fig. 1). None of these knots 
are detected in the UV images from HST/ ACS (Advance Camera 
for Surveys), acquired with F250W and F330W filters and taken 
also from the scientific archive. The wider FOV of STIS image 
(~ 28 x 51 arcsec 2 ) shows nearly 70 knots, which seem more dis- 
persed with increasing radius. Counterparts in other bands (F555W, 
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J, Ks, none detected in the L'-band) exist only for the 25 knots in 
the common FOV. In the /f-band 2 of these knots are outside of 
the NICMOS FOV. The photometry was extracted from a circular 
aperture (r a = FWHM in each filter) centred on a common position 
for all images. The background emission was then subtracted as 
the mode value for an annular concentric region (r s ~ [0'.'4-0'.'6]). 
Some of these regions stand out in the narrow-band F658N image, 
centred on Ha+[N II]. To estimate the flux in this blend, the contin- 
uum contribution was inferred by linear interpolation between the 
F555W and J broad-band filters. The possible uncertainties derived 
in the continuum estimation will be discussed later. Finally, the Ha 
flux was then corrected by the [N II] contribution by assuming that 
to be 40 per cent of the total (Baldwin et al. 1981). 



3 RESULTS 

Most of the 25 knots are unresolved or barely resolved with 
sizes of about 0'.'12 (< 10.5 pc), i.e. the spatial resolution in the 
VLT/NaCo Aj-band image. They show a median brightness of 
TflK — 20.9 mag, which corresponds to an absolute magnitude 
of Mk = —10.4 mag, too bright for being individual red super- 
giant stars (Tabur et al. 2009). The knots are distributed in the 
— 12 ^ Mk ^ —9 mag range, except knot #1, which is by far the 
brightest with Mk #1 = —13.6 mag. In the optical, they are char- 
acterized by a median mv = 24.3 mag and mi = 22.1 mag, cor- 
responding to My = —7.0 mag and Mj = —9.2 mag. Most of the 
25 knots are strong emitters in Ha. Moreover, the Ha equivalent 
width (EWHa) distribution shows > 50 A for most of them, we 
therefore adopted this value as a threshold to separate Ha emit- 
ting knots from those without Ha emission. Assuming these to 
be star clusters, their ages were derived from their EWHa. STAR- 
BURST99 models (Leitherer et al. 1999) were used, assuming an 
instantaneous burst model with the highest available metallicity, 
Z = 0.04. This follows the results by Pierce et al. (2005) who found 
higher than solar metallicity abundances (Z > 0.08) in the central 
~ lOarcsec (~ lkpc) region of NGC 1052 from the analysis of 
spectral indices measured at various radii in an optical long-slit 
spectrum. Notice that these are stellar metallicities and not nebu- 
lar ones. A lower limit for the metallicity will be considered below. 

A total of 1 5 knots were found with EWHa ^ 50 A (blue circles 
in Fig. 1), which yields ages o/<6.8Myr. Their median EWHa 
value is 116 A (6.2 Myr), with half of them located in the 105- 
175 A range (5.4-6.5 Myr), i.e. the first and third quartiles, re- 
spectively. Lower EWHa for the remaining knots may be due to 
the clusters being older or to a higher extinction which will affect 
the Ha continuum interpolation. In the conservative side, choos- 
ing larger photometric aperture in the F555W and ./-bands -the 
ones used to define the Ha continuum flux- decreases the num- 
ber of knots over the 50 A limit to 8, instead of 15. The net effect 
of considering lower metallicities, e.g. Z = 0.001, would translate 
in older ages for a given EWHa, but values over the 50 A threshold 
still limit the age to ^ 12.2 Myr. The number of "young" regions 
just increases to 16 by assuming the [N II] contribution to be 20 per 
cent instead of 40 in the Ha+[N II] blend. 

The Ha luminosity (Lua) for the 15 knots with 
EWHa > 50 A spans the [1-33] x 10 36 ergs _1 range. Notwith- 
standing, ionizing photons may be produced by an old stellar 
population due to the "UV-upturn" (see O'Connell 1999, and 
references therein). This additional contribution in the UV range 
originates in some low-mass, helium-burning stars in the extreme 
horizontal branch (EHB) and subsequent phases of evolution. In 



this scenario, the compact knots in NGC 1052 could be either 
old globular clusters (GCs) with an EHB component or young 
stellar clusters probably linked with the merger event occurred 
in the recent history of this galaxy. A detailed analysis was then 
performed for knot #6 (Fig. 1), which is representative of the 15 
knots' sample in terms of brightness and SED's shape. From its 
Ha luminosity (5.2 x 10 36 erg s _1 ) we derived an ionizing photon 
rate of ~3.8 x 10 48 pfis _1 (Osterbrock 1989, following). This 
estimate is compared in the next sections with predictions of a 
SED from an old and a young stellar population, respectively. 



3.1 Old globular clusters? 

Fig. 2 shows the spectral energy distribution (SED) of the knot #6, 
which is derived from broad-band photometry (blue dots) including 
the upper limit in the L-band. First, we show the SED of a 15 Gyr 
old, Z = 0.05 population from Bruzual & Chariot (2003, red line in 
Fig. 2) scaled to the NIR SED of knot #6 (the NIR SED comprises 
3 points: //S77F160W, NaCo J and A"s-bands). We choose scaling 
to the NIR to minimize dust extinction effects. The rate of ionizing 
photons produced by this population -mainly by blue horizontal 
branch stars- is then determined by integrating the template below 
the Lyman limit (shaded area below the red line). This results in 

1.2 x 10 47 phs -1 , i.e. more than an order of magnitude lower than 
the value measured for knot #6; thus, the Ha emission cannot be 
explained by this 15 Gyr old population. 

However, the EHB component is not included in the 15 Gyr 
old template. EHB stars present small envelopes around the helium 
core, showing a hot thermal spectrum with T > 20000 K. The in- 
tegrated energy distribution of such an old population (> 10 Gyr) 
increases shortwards of ~ 2000 A and is thus capable of emitting 
FUV photons (Greggio & Renzini 1990). To estimate the possi- 
ble contribution of these stars we consider the spectrum of the 
UV-strong elliptical galaxy NGC 4552, which is expected to have 
an important contribution from EHB stars. For consistency, this 
spectrum was scaled to the 15 Gyr old template used before in 
the optical region (yellow line in Fig 2). Shortwards of 900 A, 
NGC 4552 spectrum was extended with a T = 35000 K black body 
that was scaled to NGC 4552 in the UV part (black-dashed line in 
Fig. 2). In this case, the integrated number of ionizing photons is 
1.8 x 10 47 phs -1 , still an order of magnitude below the measured 
value. It should be noticed that the the measured Ha emission in 
the knots, and thus the ionizing photon rate, may still be an lower 
limit due to the possible contribution of Ha in absorption. This is 
expected from intermediate-age A-stars population in the clusters. 



3.2 Young stellar clusters? 

We now consider a young 6 Myr old, Z = 0.04 star cluster model 
from STARBURST99 (blue-solid line in Fig. 2, Leitherer et al. 
1999). The model choice is based on the EWHa of knot #6. Scal- 
ing it to the NIR SED region of this knot, the resulting photon 
rate is 2.5 x 10 49 phs -1 . This is more than an order of magni- 
tude larger than the measured value. Moreover, the model largely 
overestimates the emission in the optical range (see Fig. 2). As 
dust is expected in the environment of young clusters, different 
extinction values were applied to this model till matching the 
observed both SED and Ha luminosity. We started with an ini- 
tial value of Av = 2.5 mag -inferred from the Ha ratio between 
model and observation- and apply it consecutively to the SED. 
The process of fitting both Ha emission and SED shape was 
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Figure 2. Knot #6 SED (blue circles) compared with a 15 Gyr, Z = 0.05 
template (Bruzual & Chariot 2003, red line), a composite spectrum of 
NGC 4552 (Yi et al. 1998, yellow line) plus a blackbody spectrum 
(T = 35000 K, black-dashed line) and a 6 Myr, Z = 0.04 template from 
STARBURST99 (Leitherer et al. 1999, blue-solid line), also shown with an 
obscuration of Ay =3.4 mag (blue-dotted line). Shaded areas correspond 
to photons with A < 910 A. 



done in an iterative manner (using the extinction curve by Cardelli 
et al. 1989). The best match was found for a model with Ay = 
3.4 mag (blue-dotted line in Fig. 2), which reproduces the shape 
of the measured SED. The intrinsic photon rate inferred from the 
model is 4.9 x 10 49 phs~ 1 , Lhc«~6.6x 10 37 ergs _1 and mass 
1.3 x 10 4 Mq (from Ks-band luminosity, assuming an age-M/ Lk 
dependence from Leitherer et al. 1999). Note that, as we are assum- 
ing a foreground dust screen, these values may be lower limits. 



3.3 Colours 

All the 25 knots look different in brightness and colours with 
respect to those of the GCs studied in the outer 2-20 kpc of 
this galaxy by Pierce et al. (2005). The latter present a me- 
dian mv ~ 22.0 mag with V — J~ 0.9 mag, in contrast with 
the V — 7 = 2.4mag ! shown by the former (Fig. 3). Note the 
#S77F28X50LP filter is even bluer than the /-band. The differ- 
ences suggest that both populations are intrinsically distinct. This 
filter also includes the Ha+[N II] blend (~ 0.660 urn at z = 0.005) 
very close to the transmission peak 4 , but this contribution to the 
broad-band magnitude (~0.3dex for knot #6) is not enough to 
explain the large difference in V — I colour. Only knot #1 has a 
mv = 21.9 mag, compatible with that of the GCs, but its V — I 
colour is much redder (2.6 mag). 

This discrepancy between the two populations cannot be as- 
cribed to the colour bimodality reported by Forbes et al. (2001), 
since all the GCs in this study lie in the range 1.2 <V — 
I < 1.7 mag. On the other hand, an obscuration of Av = 3.4 mag 
may explain why the knots in the central region look much red- 
der than those of the GCs. Similar extinction values and colours to 
those shown by the knots have been reported for young stellar clus- 
ters (YSCs) in the starbursts galaxies NGC 5253 (Vanzi & Sauvage 
2004) and NGC 253 (Fernandez-Ontiveros et al. 2009). 



Strictly speaking, F555W-F28X50LP colour. 

http : / / www .stsci.edu/hst/stis/ design /filters/ 
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Figure 3. In grey, V — I colours for the GCs in the outer 2-20 kpc of 
NGC 1052 (Pierce et al. 2005). In blue, the F555W-F28X50LP colour for 
the 25 knots identified in the central ~2 X 2 kpc 2 . 



4 DISCUSSION 

Recent star formation is found in the centre (2x2 kpc 2 ) of 
NGC 1052 in the form of young (~ 6 Myr) and compact (~ 14 pc) 
stellar clusters with masses of ~ 10 4 Mq. The age estimation do 
not depends strongly on the stellar population model assumed. The 
total star formation rate (SFR) in the centre, derived from the Ha 
luminosity of the 25 knots, is SFR ~7x 10~ 4 M Q yr -1 (follow- 
ing Kennicutt 1998). Assuming that these clusters are spreaded 
with similar density over the whole galaxy results in an up- 
per limit for the SFR of ~ 0.01 Mq yr -1 . The resulting spe- 
cific SFR places the galaxy in the "red and death" sequence 
(SFR/M* < 10~ 12 yr" 1 , Fontana et al. 2009). Even if considering 
the extrapolation for the whole galaxy, the total mass forming these 
clusters would be less than 5 x 10 -3 per cent of the galaxy stel- 
lar mass. This suggest that the current RSF is not dominant in the 
present evolution of NGC 1052. 

Several observational signatures indicate that NGC 1052 ex- 
perienced a merger in the last 1 Gyr (c.f. Introduction). The origin 
of the knots may be associated with left over activity from that 
event. Hydrodynamical ,/V-body merger simulations (mass ratio 1 :6 
to 1:10) by Peirani et al. (2010) suggest that star formation occurs in 
several bursts after the merger event rather than being a continuous 
process (see also Kaviraj et al. 2009). This is in agreement with the 
observations: we find two very young clusters (EWHa ~ 230 A; 
~ 5.1 Myr), 4 young ones (140 EWHa < 200 A; 5.2-5.6 Myr), 
7 slightly older ones (90 < EWHa < 140 A; 5.6-6.6 Myr) and 12 
with EWHa < 90 A (> 6.6 Myr). This means that we are proba- 
bly looking at either the final stage or a minimum in the star for- 
mation activity in the centre induced by a merger event. In their 
simulations, Kaviraj et al. (2009) explore extinction values in the 
E(B — V) = — 0.5 mag range, showing that observations favour 
those models with higher extinction (E(B — V) > 0.3 mag). For 
comparison, the colour excess of knot #6 suggest higher values, i.e. 
E(B — V) ~ 1.1 mag (Av ~ 3.4 mag). 

Regarding at their spatial distribution, the knots with 
EWHa 50 A are located randomly and do not seem to follow any 
special pattern except for a weak preference for the NW-SE direc- 
tion. Moreover, we do not find any clear connection with either 
the kpc-scale radio lobes at 1.4 GHz (Fig. 1) or with the extended 
emission in X-rays detected by Chandra (Kadler et al. 2004). This 
indicates that the recent star-formation episode in the centre is not, 
at least directly, driven by nuclear activity. 

The lack of high spatial resolution CO maps does not allow us 
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to connect the knots with density enhancements in the cold mate- 
rial, but the dust lanes in the N-S direction seen in the STIS im- 
age are not apparently related with the knots. However, this scat- 
tered star formation is in line with the "widespread" star-formation 
mode found by Shapiro et al. (2010) for early-type galaxies with 
multiple and mismatched kinematic components. This is the case 
of NGC 1052, which shows 2 counter- rotating gaseous compo- 
nents decoupled from the stellar kinematics (Plana & Boulesteix 
1996). Furthermore, Shapiro et al. (2010) find that those objects 
with very low SFRs (< 0.06 Mg yr" 1 ) exhibit a more centralized 
star-formation activity, which is also favoured by recent simulations 
(Peirani et al. 2010). This might explain why Forbes et al. (2001) 
and Pierce et al. (2005) did not find young clusters in their study, 
as they covered the outskirts of NGC 1052 (2-20 kpc). Neverthe- 
less, an Ha image with a wider FOV is needed to check the young 
cluster distribution on scales larger than the inner few kpc. 

Since we determine the age of the knots on the basis of 
Ha emission, we are only sensitive to very recent star-formation 
episodes (< 10 Myr) and cannot provide a reliable age for those 
knots with EWHa < 50 A (10 out of the 25). However, these show 
very similar properties (i.e. sizes, magnitudes, colours) to those 
with EWHa > 50 A, are also randomly distributed and their SEDs 
are qualitatively similar (see Fig. 1). Results from stellar population 
models on the integrated spectrum of the central ~ 1 kpc indicate 
luminosity-weighted ages of ~ 2 Gyr (Pierce et al. 2005). How- 
ever, the contribution of the 25 knots reported in this work is not 
sufficient to "rejuvenate" the integrated emission of the central re- 
gion. In line with the scenario proposed by Shapiro et al. (2010), a 
middle-age population (10 Myr < t < 2 Gyr) is still needed, that 
would be probably related with a more intense star-formation 
episode in the past. A spatially resolved spectroscopic study of the 
central kpc region would be enlightening. 



5 FINAL REMARKS 

We propose that the knots found in the centre of NGC 1052 are 
young stellar clusters (YSCs) formed in a very recent star forma- 
tion episode, probably related with the merger event that occurred 
~ 1 Gyr ago (van Gorkom et al. 1986). This is the first detection of 
resolved YSCs in the centre of an elliptical galaxy. Their presence 
has three major implications: i) we are looking at the final stage of 
the RSF occurring in the central parsecs of NGC 1052, ii) the RSF 
takes place in the form of compact and young stellar clusters sim- 
ilar to those found in starburst galaxies, and Hi) the reddening pre- 
dicted by young stellar cluster templates [E(B — V) ~ 1.1 mag] 
should be explored in future simulations of these systems, which 
used to assume a lower value as an upper limit. 

The case of NGC 1052 is probably one of the best examples 
for the study of the recent evolution of nearby early-type galaxies 
via merging with gas-rich satellites. However, high spatial resolu- 
tion observations in dynamically younger remnants, for which we 
expect a larger number of YSCs, are desired to understand in detail 
the role of mergers in the present evolution of elliptical galaxies. 



ACKNOWLEDGMENTS 

The authors acknowledge A. Marin-Franch for his useful hints on globu- 
lar cluster physics and M. Orienti for reviewing the radio data. We also 
thank the referee for detailed comments which helped to improve the orig- 
inal manuscript. This work is partially funded by the Spanish MEC project 
AYA2007-60235. 



REFERENCES 

Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 93, 5 

Bruzual G., Chariot S„ 2003, MNRAS, 344, 1000 

Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 245 

Carlson M. N., Holtzman J. A.. Grillmair C. J., Mould J. R., Griffiths 

R. E., Ballester G. E., Burrows C. J., Clarke J. T., Crisp D., Evans R. W., 

Gallagher III J. S., Hester J. J., Hoessel J. G., Scowen P. A., Stapelfeldt 

K. R., 1999, AJ, 117, 1700 
Cole S., Norberg P., Baugh C. M., Frenk C. S., Bland-Hawthorn J., Bridges 

X, Cannon R., Colless M., Collins C, Couch W., et al. 2001, MNRAS, 

326, 255 

Cooper N. J., Lister M. L., Kochanczyk M. D., 2007, ApJS, 171, 376 
Femandez-Ontiveros J. A., Prieto M. A., Acosta-Pulido J. A., 2009, MN- 
RAS, 392, L16 

Fontana A., Santini P., Grazian A., Pentericci L., Fiore F, Castellano M., 

Giallongo E., Menci N., Salimbeni S., Cristiani S., Nonino M., Vanzella 

E., 2009, A&A, 501, 15 
Forbes D. A., Georgakakis A. E., Brodie J. P., 2001, MNRAS, 325, 1431 
Forbes D. A., Sparks W. B., Macchetto F. D., 1990, NASA Conference 

Publication, 3098, 431 
Goudfrooij P., Mack J., Kissler-Patig M., Meylan G, Minniti D., 2001, 

MNRAS, 322, 643 
Greggio L., Renzini A., 1990, ApJ, 364, 35 
Heckman T. M., 1980, A&A, 87, 152 

Jensen J. B., Tonry J. L., Barris B. J., Thompson R. I., Liu M. C, Rieke 
M. J., Ajhar E. A.. Blakeslee J. P., 2003, ApJ, 583, 712 

Kadler M., Kerp J., Ros E., Falcke H., Pogge R. W., Zensus J. A., 2004, 
A&A, 420, 467 

Kaviraj S., 2010, MNRAS, pp 1161 — h 

Kaviraj S., Peirani S., Khochfar S., Silk J., Kay S., 2009, MNRAS, 394, 
1713 

Kaviraj S., Schawinski K., Devriendt J. E. G., Ferreras I., Khochfar S., 
Yoon S., Yi S. K., Deharveng J., Boselli A., Barlow T., et al. 2007, ApJS, 
173,619 

Kennicutt Jr. R. C, 1998, ARA&A, 36, 189 

Leitherer C, Schaerer D., Goldader J. D., Gonzalez DelgadoR. M., Robert 
C, Kune D. F, de Mello D. F, Devost D., Heckman T. M., 1999, ApJS, 
123, 3 

Lopez-Sanjuan C, et al., 2010, arXiv: 1009.5921 
O'Connell R. W., 1999, ARA&A, 37, 603 

Osterbrock D. E., 1989, Astrophysics of gaseous nebulae and active galac- 
tic nuclei. University Science Books 

Peirani S., Crockett R. M., Geen S., Khochfar S., Kaviraj S., Silk J., 2010, 
MNRAS, 405, 2327 

Pierce M., Brodie J. P., Forbes D. A., Beasley M. A., Proctor R., Strader 
J.. 2005, MNRAS, 358, 419 

Plana H., Boulesteix J., 1996, A&A, 307, 391 

Pogge R. W., Maoz D., Ho L. C, Eracleous M., 2000, ApJ, 532, 323 
Rampazzo R., Marino A., Tantalo R., Bettoni D., Buson L. M., Chiosi C, 

Galletta G., Griitzbauch R., Rich R. M., 2007, MNRAS, 381, 245 
Shapiro K. L., et al., 2010, MNRAS, 402, 2140 
Sofue Y., 1993, PASP, 105, 308 
Stetson P. B., 1987, PASP, 99, 191 

Strateva I., Ivezic Z, Knapp G. R., Narayanan V. K, Strauss M. A., Gunn 
J. E., Lupton R. H., Schlegel D., Bahcall N. A., Brinkmann J., et al. 2001, 
AJ, 122, 1861 

Tabur V., Kiss L. L., Bedding T. R., 2009, ApJ, 703, L72 

van Gorkom J. H., Knapp G. R., Ekers R. D., Ekers D. D., Laing R. A., 

Polk K. S., 1989, AJ, 97, 708 
van Gorkom J. H., Knapp G. R., Raimond E., Faber S. M., Gallagher J. S., 

1986, AJ, 91,791 
Vanzi L., Sauvage M., 2004, A&A, 415, 509 

Williams B. F., Dalcanton J. J., Gilbert K. M., Stilp A., Dolphin A., Seth 

A. C, Weisz D., Skillman E., 2010, ApJ, 716, 71 
Yi S., Demarque P., Oemler A. J., 1998, ApJ, 492, 480 



6 Ferndndez-Ontiveros et al. 

This paper has been typeset from a Tj5£/ ETgX file prepared by the 
author. 



